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ABSTRACT 

Proton acceleration at a parallel coronal shock is modeled with self-consistent 
Alfven wave excitation and shock transmission. 18 - 50 keV seed protons at 0.1% 
of plasma proton density are accelerated in 10 minutes to a power-law intensity 
spectrum rolling over at 300 MeV by a 2500km s -1 shock traveling outward from 
3.5r 0 , for typical coronal conditions and low ambient wave intensities. Interac- 
tion of high-energy protons of large pitch-angles with Alfven waves amplified by 
low-energy protons of small pitch angles is key to rapid acceleration. Shock ac- 
celeration is not significantly retarded by sunward streaming protons interacting 
with downstream waves. There is no significant second-order Fermi acceleration. 

Subject headings: acceleration of particles — shock waves — Sun: coronal mass 
ejections (CMEs) 


1. INTRODUCTION 

The acceleration of solar energetic particles (SEPs) by coronal mass ejection (CME) 
driven shocks is an outstanding problem in space weather and astrophysics. Observations 
indicate that SEPs are accelerated to GeVs in < 10 minutes after CME launch in many 
events. This short time scale requires mean free paths A < 1 x 10 -6 AU in diffusive shock 
acceleration models, in apparent contradiction with A > 0.2 AU deduced from interplane- 
tary (IP) transport of SEPs (e.g., Bieber et al. 1994). Waves excited by streaming SEPs 
resolve this paradox. Self-amplified Alfven waves were invoked in the confinement (Wentzel 
1974) and diffusive shock acceleration of cosmic rays (Bell 1978). Steady-state models with 
self-consistent Alfven waves were applied successfully at the Earth’s bow shock and traveling 
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IP shocks (Lee 1982, 1983, 2005; Gordon et al. 1999). Ng & Reames (1994) showed theoret- 
ically that self-amplibed waves limit SEP intensity at 1 AU in large events. Reames & Ng 
(1998) inferred streaming-limited SEP intensities and their radial dependence from Helios 
and GOES data. Ng et al. (1999) and Reames et al. (2000) showed how proton-amplified 
waves control abundance variation in SEP events. Can streaming SEPs amplify Alfven waves 
upstream of CME-driven coronal shocks by a factor > 10 5 in a few minutes and to what 
energy are the SEPs accelerated? The answer requires a time-dependent model that treats 
sclf-consistently the nonlinear, coupled evolution of SEPs and Alfven waves. 

SEP acceleration at a CME-driven shock is a very complex phenomenon (Lee 2005; 
Cliver et al. 2004; Reames 1999). Researchers have studied different aspects theoretically 
with various simplifying assumptions, using various analytical and numerical techniques. 
Significant progress has been made with many issues under lively debate. The reader is 
referred to Tylka & Lee (2006), Zank et al. (2006), and Giacalone (2005) for the effects of 
shock geometry; Zank et al. (2007) for the excellent work performed by their group; and 
Giacalone & Neugebauer (2008) for the effects of a rippled shock. 

In this Letter we address the question of how rapid and to what energy a fast parallel 
shock accelerates SEPs in typical coronal conditions by solving numerically the coupled 
evolution of SEPs and Alfven waves. We also address the issues of whether shock acceleration 
is significantly retarded by SEP damping of outward waves and excitation of inward waves 
downstream of the shock (Vainio 2001; Ng 2007) and whether second-order Fermi acceleration 
plays a significant role. 

Our model includes in SEP transport: focusing, convection, adiabatic deceleration, 
pitch-angle scattering, and momentum diffusion due to scattering by counter-streaming 
Alfven waves. Treatment of Alfven waves includes propagation, SEP-driven growth and 
damping, and shock transmission (Vainio & Schlickeiser 1999). A unique feature of the 
numerical model is the implementation of the pitch-angle dependent cyclotron resonance 
condition to calculate self-consistent SEP scattering and Alfven wave growth. This feature 
is essential for understanding how acceleration is initiated for SEPs above the “knee” energy 
Eknee ■ Momentum diffusion governed by the diffusion tensor elements D^, D^p = Dp 
Dpp (Ng et al. 2003) is solved by using an iterative method more accurate than the direct 
method employed in Ng (2007). The model now covers larger rigidity (P) and wavenumber 
( k ) ranges to reach higher energy. 

Our model accounts for wave evolution time and does not adopt an instantaneous steady- 
state upstream (US) wave spectrum (Li et al. 2005). The model of Vainio & Laitinen (2007) 
treats US wave growth and accelerates protons to tens of MeV in 10 minutes. They exclude 
the downstream (DS) region by reflecting SEPs from the shock statistically. The DS region 
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is an important part of onr consideration. 

We report results for a strong 2500 km/s shock starting at 3.5r 0 . Oliver et al. (2004) have 
reported that shock acceleration is most efficient above ~ 3r 0 . Studies on the dependence 
on shock, plasma, and turbulence parameters will be published in a full paper. We describe 
the model in Section 2, present the results in Section 3, and close with a discussion in Section 
4. 


2. THE MODEL 


The time-dependent nonlinear model of Ng et al. (2003) is generalized to treat SEP 
shock acceleration and shock transmission of Alfven waves in a computational box comoving 
with the shock. Radial geometry is assumed with mean magnetic held B = B 0 (r 0 /r) 2 , 
plasma proton density n H = n^ 0 (ro/r) 2 , hence Alfven velocity V A = V Afi (r 0 /r), where B 0 , 
n H , o, and U A , 0 = B 0 / ^/47m H , 0 are values at the reference radius r 0 = 3.5r 0 . The evolution of 
the energetic proton distribution is governed by: 
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where F(P, /i, x,t) = ( B 0 /B)(P/P 0 ) 3 f(P,n,x,t ); / = proton phase-space density in mixed 
coordinates: distance x from shock and time t in the Sun’s inertial frame; P, proton velocity 
v, and pitch-angle cosine yU in the local outward wave frame moving with velocity W + {x,t) = 
Wv + V A relative to the Sun. dW + /dt = <9*14^ +xd x W + ] y = (1 + /itTly/c 2 ) -1 ; q = ln(P/P 0 ); 
P 0 = constant; D /ig = D m = P^D^p, D ee = P^ 2 P pp ; r(x,t) = x + r sh (t); r sh (t) = shock 
radius; E s h = constant shock velocity. The plasma velocity V^rr) = V ( j at x < 0 (DS) and 
V sw (x) — V u at x > 0 (US). Besides focusing and adiabatic deceleration, equations (3) and 
(4) include the differential frame transformation of (/x, P) following particle motion. The 
shock jump in W^(r,t) gives a P-increment for almost all shock-crossing particles (eq. [4]). 


We consider only right- and left-hand circularly polarized, inward and outward parallel 
propagating Alfven waves. The ambient outward wave magnetic intensities / R+ = / L+ oc k~ s 
are specified for US and DS separately via steady-state solutions to the wave kinetic equation 
(Ng et al. 2003). As outward waves dominate near the Sun, we specify J R " = J L ~ = 0.05/ R+ . 
I a , n H and V A are then modified for the shocked DS region assuming 10 s shock transmission 
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before seed-particle injection at t — 0. Wave evolution is governed by the wave kinetic 
equation 

d t * ff + d x [(W ff - V sh )^] + vdr,(Vf* ff ) = 7^, (5) 

with wave action density ^ a (r},x,t) = 2I a (r},x,t)rj\W <J \/(r}oV A ), g = k/B , r/ 0 = constant, 
17 CT = 2W a /r — dW a /dr, = fractional growth rate of / CT , HA = W + for a = P+, L+, and 
W a = HA_ = V sw (r,t) — V A (r,t ) for cr = R—,L—. Equations (1) and (5) are coupled via 

fdkR^r, (6) 

cr cr 

7<r(fc, r, t) = 2n 2 ce 3 g a V A Jj djulPR^Gf, (7) 

where Gf = P 3 [£(l - /il4/u)]^ 2 [(9 /t / - {V a /v){iid^f - d e f)\; £ = total particle energy; 
g a = 1, V a = 0 for cr = P+,L+; and g a = —1, V a = — 2V A for a = R—,L~. The wave- 
particle resonance function R° (fj,, P, k, V A , B) (Ng & Reames 1995; Ng et al. 2003) includes 
minimal broadening of the quasilincar resonance condition k = B/[P(/j, — V a /v)\. 

At the initial shock location at r = 3.5 r Q , we specify V u = 83 km s _1 , n H , 0 = 2 x 10 5 
cm -3 , V A 0 = 700 km s _1 from the semi-empirical models of Guhathakurta et al. (1999). 
Hence B 0 = 0.143 gauss, consistent with results from gyrosynchrotron emission (Dulk & 
McLean 1978). Prescribed ambient I a oc /c _L5 gives A = 0.23 AU at 1 MeV (0.68 AU at 
100 MeV). Constant V u , V ( \ = 1880 km s -1 , and Vh = 2500 km s -1 in the Sun’s inertial 
frame are assumed, giving a fluid compression ratio c/ = 3.895 and an initial Alfven Mach 
number Ma = 3.45 increasing to 5.56 at t — 10 min, when r s h = 5.65 r Q . The relation 
between cy, M Al plasma beta = 0.1, and relative wave amplitude (Vainio & Schlickeiser 

1999) is satisfied at t — 0 but not at t > 0. The back reaction of SEPs and amplified waves 

on shock evolution is not considered. 

For quasi-parallel shocks, Zank et al. (2001) calculate that 2% of the solar-wind protons 
are reflected by cross-shock electric potential to > 3P ramp . Here 3P ramp = 90 keV. From the 
observed post-shock proton distribution in the August 27, 1978 event (Gosling et al. 1981, 
Figs. 4 and 6), we find tail density above 18 keV is a fraction b = 0.0002 of the observed 
n H = 25 cm -3 . Larger b is expected in the strongly heated plasma DS of the 2500 km s^ 1 
shock. There may also be suprathcrmal remnants from previous SEP events (Mason et al. 
1999; Tylka et al. 2001). Guided by the above, we prescribe DS isotropic seed protons 

/»J = )(a - 3)f>-«/{47r(Py« - (8) 

at P a < P < Pb, with injection fraction b = 0.001, a = 5, P a = 6 MV (18 keV), and 
P b = 9.8 MV (50 keV). Note that P b < P\ = 10.4 MV, the lowest grid rigidity, so the DS 
seeds emerge US as outward proton beams from 10.4 to 14.7 MV (57 to 115 keV). 
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Each US wave I u {k u ) crossing the shock is converted to two DS waves Id{kd ) of the same 
hclicity, taking account of shock amplification and wavenumber increase: 

I d (k d )/I u (k u ) = (1 + ( u ( d /y/cj) 2 c w /4: (9) 

with wave compression ratio c w = kd/k u = c f{M a~ Cu) / {M a~ Cd^/c]) , Cu — ±1 and Q = ±1 
for outward (inward) US and DS waves (Vainio & Schlickeiser 1999; Webb et ah 1999). Since 
wave growth is limited by nonlinear processes, we impose I a {k) < I sa t(k) = B 2 /{‘ink) so that 
A > 3 r g . Shock transmission of US outward waves produces only a small fraction (5%) of 
DS inward wave intensity. In applying the above limit, we preserve the ratio of transmitted 
inward and outward waves. 

The numerical grid is: Xj = 2.15 x 10 _4 (j — 1/2) r Q , j = -200:800; ji t = 0.05 (i+1/2), 
i— -20:19; Qe = Qi + 0.08664(1' — 1), t = 1:60, P\ = 10.38MV; and log rj m = log 771 + 
0.05(m — 1), m = 1:58, rj\ = 5.465 x 10~ 3 4 MV _1 . At the boundaries, exterior F = interior 
F. At the outer boundary, exterior I a = interior I a . The /x and P increments of shock- 
crossing protons are evaluated by frame transformation. Equations (1) and (5) are solved 
with dynamic time steps using operator splitting and locally one- dimensional differencing 
(Ng et al. 2003); (/x, P)-diffusion is solved iteratively. 


3. RESULTS 

In 600 s, the 2500 km s _1 parallel shock accelerates 18 - 50 keV DS seed protons to 
a power-law intensity spectrum j E extending to E knee ~ 300 MeV (Fig. la). The energy 
spectral index of 0.8 is consistent with c w = 6.46 decreasing to 4.95 between US and DS 
outward Alfven waves. At 12.3 MeV, the j E spatial profile rises swiftly from a spike at 
180 s to the steady-state form (Lee 1983): a DS plateau and an US exponential tail at 
t > 250 s (Fig. lb). This implies negligible second-order Fermi acceleration, consistent with 
no enhanced inward and outward waves in proximity. 

Streaming SEPs US rapidly amplify outward waves and damp inward waves. Wave 
growth extends to low wavenumber k (Fig. 2a) as acceleration proceeds to high energy. The 
amplified waves enhance scattering, providing a positive feedback that rapidly bootstraps 
acceleration. The DS outward waves derive from shock-transmitted US outward waves and 
are shifted to higher k (Fig. 2b; cf. Fig. 2a). The much weaker DS inward waves (not 
shown) are h-shifted by a smaller amount. The h-shifted enhanced DS waves thus resonate 
with SEPs of lower values of P{fi — V a /v). The effect is seen by comparing D ^ US and DS 
at t = 576 s at five rigidities from 140 to 1449 MV (Figs. 2c, 2d). The hump in each D 
curve originates from proton-amplified US waves and occupy a smaller /x range DS than US 
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at P > 1117 MV. At P < 305 MV, however, the DS waves at k resonant with large // are 
amplified by US SEPs of higher P, so the DS humps are also wide. 

How acceleration begins can be seen in the evolution of US and DS //-distributions 
(Figs. 2e, 2f). The scale factors chart the course. SEPs at E > Eh nee are first accelerated at 
the D IHl hump at small |//|, where the resonant waves are available. Some accelerated SEPs 
US scatter to larger // and excite US waves at smaller k, which are then transmitted DS. 
The result is that the DS D hump widens and acceleration extends to higher energy. The 
//-dependent resonance condition k ~ B/[P(fi — V a /v)\ is crucial in this process. It allows 
waves excited by low P large // SEPs to scatter high P small // SEPs (and vice versa). At 
t — 0, I H+ (k), I h+ (k) spectra are already shock-transmission enhanced DS relative to US by 
a factor of 60 (Figs. 2a, b). Streaming SEPs do not always amplify (damp) waves traveling 
in the same (opposite) direction, because the //-distribution is not always monotonic (Figs. 
2e,f). Wave growth and damping do eventually extend to lower k. 

Do DS inward waves, which scatter SEPs back to the shock with energy loss as seen 
in outward-wave frame, dominate in the immediate DS and retard acceleration? Figure lc 
compares D+ and D~ due to outward and inward waves, respectively, of 48.2 MeV protons 
in the first DS cell at t — 72 and 576 s. At both times, scattering by strong outward 
waves dominate in the central //-range, where most of the SEPs are scattered back and forth 
across the shock. Hence shock acceleration is not significantly retarded by SEP damping 
(amplification) of DS outward (inward) waves. 


4. SUMMARY AND DISCUSSION 

Shock acceleration is simulated in a moving box model with the usual processes of SEP 
transport and momentum diffusion, and the propagation, shock transmission, and SEP- 
driven evolution of Alfven waves. For typical coronal conditions and weak ambient waves, 
a 2500km s -1 parallel shock starting at 3.5 r 0 accelerates 18-50 keV DS seed protons at 
0.001n H to a power-law intensity spectrum with E^nee ~ 300 MeV in 10 min (Fig. la). The 
box size of O.215r 0 = 2.8 x 10 5 ion inertial lengths and fit = 8.2 x 10 5 for t = 600 s are 
respectively 10 and 140 times larger than in the hybrid simulation of Giacalone (2004). 

A key feature of the model is the consistent use of the //-dependent cyclotron resonance 
condition B/k ~ P(/z — V A /v) to calculate the momentum diffusion tensor D a p and the wave 
growth rate 7 CT , so that waves excited by low-energy protons also scatter high-energy ones. 
Using the //-dependent condition to calculate 7 CT but the “sharpened” condition B/k = P to 
calculate D underestimates the effect of amplified waves on D^, giving little acceleration 
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(Berezhko et al. 1998). Using the “sharpened” condition to calculate both 7 ^ and D /t/i (Vainio 
& Laitinen 2007) requires high-energy protons to amplify unique resonant waves from the 
ambient state before they experience enhanced scattering. The “sharpened” condition is a 
compromise on the physics of cyclotron resonant interaction. 

The evolution of SEPs and DS inward and outward Alfven waves shows that (a) second- 
order Fermi acceleration is not significant, and (b) in the immediate DS, D+ » D~^ 
(Fig. lc) over the //-range where initial acceleration at E > Eknee occurs, hence SEP inter- 
action with DS waves do not seriously retard acceleration. 

At t — 600 s, SEP energy density e SEP = 7900 MeV cm -3 at r s h = 5.65r©, falling steeply 
with x; shock ramp energy density e ramp = 6100 MeV cm -3 , Alfven- wave magnetic energy 
density e wave = 470 MeV cm -3 , and Bg/8n = 510 MeV cm -3 . This indicates the shock 
would be modified, an effect not considered here, and that we are pushing the limits of 
quasilincar theory. Reducing injection from b = 0.001 to b = 0.0005 yields E knee ~ 100 MeV, 
e SEP = 1950 MeV cm -3 and e wave = 420 MeV cm -3 . Unfortunately, it is impossible to deduce 
Je at 1 AU from Figure la because of streaming limit (Ng & Reames 1994; Reames & Ng 
1998). We plan to inject escaping proton flux from this model into an IP transport model 
with self-consistent waves (e.g., Ng et al. 2003) to predict ;je at 1 AU. 

Detailed study of model dependence on shock, ambient wave, and plasma parameters 
will be presented separately in a full paper. Although the model is simplified, neglecting e.g. 
shock geometry and wave cascading, it captures the essential role of wave-particle interaction 
with a realistic description of SEP shock acceleration through self-amplified waves. Faster 
acceleration to higher energy may be possible at quasi-perpendicular shocks. In future, we 
will include minor ions and the effects of shock geometry, and integrate the box acceleration 
model into an IP transport model. 


This work was supported by NASA under LWS-04-0000-0076 and SHP04-0016-0024. 
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Fig. 1. — (a) Proton intensity spectra at the first cell upstream, (b) 12.3 MeV proton 
intensity spatial profiles, (c) 48.2 MeV proton /^-diffusion coefficients D± vs /j due to outward 
and inward Alfven waves at the first cell downstream at t = 72 s and 576 s. Ap = 1/20. 
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Fig. 2. — Evolution in the first cell upstream of (a) J R+ at 72s intervals, (c) D MJ at t = 576s 
at indicated rigidities, and (e) pitch-angle distributions at indicated times. The lower panels 
(b), (d), and (f) show the same respectively in the first cell downstream. 


